Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 1 February 2008 (MN MfeK style file vl.4) 



00 

On 
On 



The Angular Correlation Function of K' ~ 19.5 Galaxies 
and the Detection of a Cluster at z = 0.775 

Nathan Roche 1 ' 4 , Stephen A. Eales 1 ' 5 , Hans Hippelein 2 ' 6 
and Chris J. Willott 3 ' 7 

1 Department of Physics and Astronomy, University of Wales Cardiff, P.O. Box 913, Cardiff CF2 3YB, Wales. 
2 Max Planck Institut fur Astronomie, Konigstuhl 17, 69117 Heidelberg, Germany. 

3 Department of Astrophysics, University of Oxford, Nuclear and Astrophysics Laboratory, Keble Road, Oxford OX1 3RH, England. 
4 ndr0astro.cf.ac.uk 5 sae@astro.cf.ac.uk 6 hippeleiOmpia-hd.mpg.de 7 cwillottlOphysics.ox.ac.uk 

1 February 2008 

ABSTRACT 



We investigate (i) the clustering environment of a sample of 5 radio galaxies at 

■ 0.7 < z < 0.8 and (ii) the galaxy angular correlation function, uj(9), on five -KT'-band 

(2.1/im) images covering a total of 162.2 arcmin 2 to a completeness limit K' ~ 19.5, 

Applying two methods - counting galaxies within 1.5 arcmin of each radio galaxy, 
and using a detection routine with a modelled cluster profile - we detect a cluster 
of estimated Abell richness = 85 ± 25 (class 1 or 2), approximately centred on 

■ the radio galaxy 5C6.75 at z = 0.775. Of the other radio galaxies, two appear to be 
qq | in less rich groups or structures, and two in field environments. The mean clustering 
Q\ . environment of all 5 radio galaxies is estimated to be of Na — 29± 14 richness, similar 

to that of radio galaxies at more moderate redshifts of 0.35 < z < 0.55. 

The angular correlation function, u>(0), of the detected galaxies showed a positive 
signal out to at least ~ 20 arcsec, with a ~ 4tr detection of clustering for magnitude 
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O ' limits K' = 18.5-20.0. The relatively high amplitude of uj(6) and its shallow scaling 

with magnitude limit are most consistent with a luminosity evolution model in which 
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CZ) ■ E/SO galaxies are much more clustered than spirals (rrj = 8.4 compared to 4.2 h 

Mpc) and clustering is approximately stable (e ~ 0) to z ~ 1.5, possibly with an 
increase above the stable model in the clustering of red galaxies at the highest (z > 1.5) 
redshifts. 

rN \ Our images also show a significant excess of close (1.5-5.0 arcsec separation) pairs 

^ . of galaxies compared to the expectation from uj(9) at larger separations. We estimate 

that a 11.0 ± 3.4 per cent fraction of K' < 19.5 galaxies are in close pairs in excess 
of the observed oj(9), if this is of the form oj(9) oc 8~ Q - S . This can be explained if the 
local rate of galaxy mergers and interactions increases with redshift as ~ (1 + z) m 
with to = 1.33±gJ5- 
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1 INTRODUCTION 

Statistical measures of the clustering of distant galaxies, and 
the identification of rich galaxy clusters at high redshifts, are 
both of great importance in the study of galaxy evolution 
and the formation of structure in the Universe. The clus- 
tering of galaxies on the sky, as described in terms of the 
angular correlation function, u)(6), has been studied exten- 
sively at A ~ 0.4-0.9 /im wavelengths. The amplitude of u)(6) 
decreases on going faintward approximately as expected if 



galaxies undergo moderate luminosity evolution while their 
intrinsic clustering remains approximately stable in proper 
co-ordinates (e.g. Roche et al. 1996; Roche and Eales 1998; 
Postman et al. 1998). 

It is now possible to survey sufficiently large areas at 
near infra-red wavelengths, e.g. the 2.2/xm A"-band, to mea- 
sure the galaxy u{0) amplitude (Baugh et al. 1996; Carl- 
berg et al. 1997). Near-IR surveys, compared to those at 
visible light wavelengths, contain a much larger proportion 
of early-type galaxies and are less sensitive to the effects of 
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increased star-formation in the past, so may be particularly 
useful in studying the evolution of clustering separately from 
that of the star-formation rate. Previously (Roche, Eales 
and Hippelein 1998, hereafter Paper I), we measured the 
uj{9) amplitude to K — 20 on 17 small fields totalling 101.5 
arcmin 2 , imaged with either the Redeye near-infra-red cam- 
era on the Canada-France Hawaii Telescope or the Magic 
camera on Calar Alto. The u){9) amplitude appeared to be 
relatively high compared to that from blue-band galaxy sur- 
veys of comparable depth (e.g. Roche et al. 1996), suggesting 
the intrinsic clustering of red E/SO galaxies is significantly 
stronger (r ~ 8 hT 1 Mpc) than that of spirals and that 
clustering is approximately stable (e ~ 0) with redshift, but 
more K-bamd data was needed to improve the statistics. 

Some faint galaxy surveys have shown an excess in the 
number of close pairs of galaxies compared to the number 
expected from the inward extrapolation of ui(9) at larger 
scales (e.g. Infante, de Mello and Menanteau 1996). The 
excess pairs appear to be interacting or merging galaxies, 
and the number observed suggest some increase with red- 
shift in the merger/interaction rate, an important finding as 
the merger rate evolution is of cosmological interest (Carl- 
berg, Pritchet and Infante 1994). However, there have been 
large differences between the close pair statistics of different 
datasets. For example, in Paper I, we found some excess of 
2-3 arcsec separation pairs on the 12 fields observed in good 
seeing with the Redeye camera, but not on the data from 
the Magic camera, which was of lower resolution due to both 
poorer seeing and a larger pixelsize. Similarly, our R < 23.5 
Wide Field Camera survey (Roche and Eales 1998) found 
the u>(0) amplitude at 2 < 6 < 5 arcsec to be significantly 
higher that at larger 6, indicating an excess of close pairs, on 
one 1.01 deg 2 field but not on a second 0.75 deg 2 field, where 
again the seeing was poorer. These inconsistencies suggest 
that faint galaxy close pair counts are critically dependent 
on data quality, and require similar analyses of further deep 
surveys. 

The identification of rich clusters of galaxies at high red- 
shifts (z > 0.5) is also of great interest - their number, rich- 
ness and other properties may provide important constraints 
on cosmological parameters (e.g. Eke, Cole and Frenk 1996). 
At visible light wavelengths, the contrast of clusters against 
the field galaxies falls steeply with increasing cluster red- 
shift, making detection difficult at z > 0.6, but this is in 
part due to the strong k-correction dimming of early-type 
galaxies. In the K-band, the k-correction produces a bright- 
ening to z ~ 1, for all galaxies, so clusters should be de- 
tectable to higher redshifts. Near-IR surveys have already 
found rich clusters at redshifts as high as z — 1.27 (Stanford 
et al. 1997) and possibly at z = 2.39 (Waddington 1998). 

Searches for high redshift clusters have often concen- 
trated on areas containing high redshift QSOs or radio galax- 
ies, as some of these sources do lie in rich clusters and the 
AGN provides an easily measured redshift. The clustering 
environment of radio galaxies is of interest in itself for un- 
derstanding their evolution. At low redshifts the more lu- 
minous (FRII) radio galaxies are generally found in field 
environments, whereas at z ~ 0.5, luminous radio galaxies 
are found in similar numbers in the field and in rich clusters, 
with the average environment being approximately that of 
an Abell class cluster (Yates, Miller and Peacock 1989; 
Hill and Lilly 1991). 



Yee and Ellingson (1993) explained this as an environ- 
mental influence on radio galaxy evolution, whereby sources 
in clusters decrease much more rapidly in radio luminosity at 
z < 0.5 than those in the field. Wan and Daly (1996) found 
no observable differences between the powerful radio galax- 
ies in field and cluster environments, and concluded that 
the change in radiogalaxy environment with redshift is most 
likely due to an evolution of the clusters themselves - an in- 
crease in the typical intra-cluster medium pressure between 
z ~ 0.5 and z ~ might surpress FRII activity in most rich 
clusters at lower redshifts and allow only the less radiolumi- 
nous FRI outbursts to occur. However, to better understand 
this process, it is important to determine whether there is 
any further change in the mean radio galaxy environment 
at z > 0.5. 

In Paper I, 17 high redshift (z mean = 1.1) 6C radio 
galaxies were cross-correlated with the surrounding K < 20 
galaxies to estimate their mean clustering environment. No 
significant signal was detected in the cross-correlation func- 
tion - the upper limits were consistent with a mean radio 
galaxy environment similar to that at z ~ 0.5 (i.e. Abell 
clusters) but argued against a much richer environment (i.e. 
Abell 2 clusters). However, at similar redshifts, some of the 
3CR radio galaxies do inhabit rich clusters e.g. 3CR184 at 
z = 0.996 (Deltorn et al. 1997) and 3C336 at z = 0.927 
(Bower and Smail 1997). 

In this paper we investigate the clustering environment 
of 5 radio galaxies at 0.7 < z < 0.8, observed in the K'-band 
with the large format Omega camera on Calar Alto (see 
Section 2.1). The more moderate radiogalaxy redshifts and 
larger field areas compared to the Paper I survey should im- 
prove our chances of finding clusters. Section 2 describes the 
observational data and the detection of galaxies and stars, 
Section 3.1 the visual impression of the radio galaxy envi- 
ronments and Section 3.2 the modelling of distant clusters. 
We search for clusters centred on the radio galaxies using 
two methods, described in Sections 3.3 and 4 respectively. 
In Section 5 we investigate the galaxy u>(9), and in Section 
6 the number of close pairs of galaxies. Section 7 discusses 
the implications of all the results. 



2 DATA 

2.1 Observations 

Our dataset consists of images of 5 fields, each centred on 
a radio galaxy at 0.7 < z < 0.8 (Table 1). One source is 
from the 6C catalog of radio galaxies with 151 MHz fluxes 
of at least 2.2 Jy, the other four from the 7C catalog (Willott 
et al. 1998, and in preparation) with a fainter flux limit of 
0.5 Jy (these have names beginning '5C, as they were first 
detected in a 5C survey). These 7C galaxies have radio lumi- 
nosities L(151MHz) ~ 10 27 ' 6 WHz -1 , much lower than 3C 
galaxies but above the FRII/FRI divide. The five galaxies 
were selected from the catalogs solely on the basis of being in 
the desired redshift range, so should be an unbiased sample 
of 0.7 < z < 0.8 radio galaxies. 

The fields were observed in the 2.1/im K'-band (Wain- 
scoat and Cowie 1992) on the nights of 13 and 14 December 
1997, using the Omega wide-field near-infrared camera at 
the prime focus of the 3.5m Calar Alto telescope (Sierra de 
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Table 1. Positions (equinox 2000.0), spectroscopic redshifts and 
K' magnitudes (as measured from this data) of the five radio 
galaxies 



Galaxy 


R.A. 


Dec. 


z 


K> 


5C6.25 


02:10:24.45 


+34:10:46.0 


0.706 


17.44 


5C6.29 


02:11:05.85 


+32:56:44.2 


0.720 


16.61 


5C6.43 


02:12:02.69 


+34:02:18.1 


0.775 


17.95 


5C6.75 


02:14:01.20 


+30:26:14.3 


0.775 


17.05 


6C0822 


08:25:47.36 


+34:24:26.5 


0.700 


17.48 



los Filabres, Andalucia, Spain). This has a 1024 x 1024 pixel 
HgCdTe (1-2.5 (im) array with pixelsize 0.3961 arcsec. Each 
field was exposed for 30 seconds 120 times, and the 120 ex- 
posured stacked together and flat-fielded during the observ- 
ing run. For the purposes of flat-fielding, the 120 exposures 
were slightly offset in a grid pattern, reducing the areas of 
the final images (in which all pixels have a full 60 minutes 
exposure) to approximately 870 x 870 pixels or 5.75 x 5.75 
arcmin. 
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2.2 Source detection 

The Starlink PISA (Position, Intensity and Shape Analysis) 
package, developed by M. Irwin, P. Draper and N. Eaton, 
was used to detect and catalog the objects. The chosen de- 
tection criterion was that a source must exceed an intensity 
threshold of 1.5<r above the background noise, or approxi- 
mately 21.0 K' mag arcsec -2 , in at least 5 connected pixels. 
As in Paper I, PISA was run with deblend (objects connected 
at the detection threshold may be split if they are separate 
at a higher threshold) and the 'total magnitudes' option. 
There were a total of 3160 detections. 

Detections in four circular 'holes' (radius 9 to 20 arc- 
sec) around bright, saturated stellar images were excluded 
in order to remove spurious noise detections, leaving a to- 
tal 162.3 arcmin area of usable data . Radial profiles were 
fitted to several non-saturated stars on each image using 
'pisafit', and from these profiles the average resolution on 
the reduced images was estimated as FWHM — 1.13 + 0.06 
arcsec. Star-galaxy separation was performed using plots of 
central against total intensity, normalized to the ratio from 
the fitted stellar profile. On these plots the stellar locus was 
separable from the galaxies to K' ~ 17, and a total of 185 de- 
tections to this limit were classed as stars. All fainter objects 
were assumed to be galaxies, but we later apply corrections 
to our results for the effects of faint star contamination. 

2.3 Galaxy counts 

Figure 1 shows the differential number counts, in A(K') = 
0.5 mag bins, of objects classed as galaxies on these im- 
ages, with field-to-field error bars. These are compared with 
galaxy counts from Paper I and a number of other K-band 
surveys, plotted here assuming a mean colour K' = if +0.13, 
and non-evolving and pure luminosity evolution (PLE) mod- 
els, which assume go = 0.05 and are the same as those of 
Roche and Eales (1998) but computed in the K'-bsmd rather 
than the i?-band. 

To K' ~ 19.5, our galaxy counts agree well with those 
from previous surveys, including those using much deeper 



Figure 1. Differential galaxy number counts, in 0.5 mag intervals 
of K' magnitude, for our five fields compared to our Paper I and 
the i4"-band surveys of Gardner, Cowic and Wainscoat (1993), 
Djorgovski ct al. (1995), Moustakas ct al. (1997) and Szokoly 
ct al. (1998), together with the predictions of go = 0.05 PLE 
(dashed) and non-evolving (dotted) models. 



data, suggesting that our galaxy detection is virtually com- 
plete. Our counts level out at 19.5 < K' < 20.0, where on 
the basis of the count gradient from models and deeper data 
( ^flgg^Q ~ 0.30) they are ~ 32 per cent incomplete, and fall 
off at K' > 20. The K and if'-band number counts are 
generally in good agreement with the PLE model over the 
entire ~ 10 mag range. 

Figure 2 shows the PLE model for the redshift distri- 
bution N(z) at the K' = 19.5 limit. The predicted mean 
redshift is 0.94, the radio galaxy redshifts are close to the 
peak of N(z), and we expect to see galaxies to a maximum 
redshift z ~ 2.4, with E/SOs predominating at z > 1. 



2.4 Star counts 

Figure 3 shows the number counts of objects classed as stars, 
with field-to-field errors. The field-to-field scatter in star- 
counts is only ~ y/N as the 5 fields are at similar distances 

o 

from the Galactic plane, \b\ = 25.9-33.5 . 

To estimate the contamination of our galaxy sample by 
the unclassified K > 17 stars, we fit the observed star counts 
with a model and extrapolate faintwards. To K' ~ 16.5, the 
slope of the star counts is close to ^ — 0.2, but K- 
band star-count models (Glazebrook et al. 1994; Minezaki 
et al. 1998) flatten to d(l °^ N) ~ 0.1 at 17 < K < 20. The 
model of Minezaki et al. (1998) is for the Galactic pole and 
gives too shallow a count at 14 < K < 17 to fit our data, 
and the Glazebrook et al. (1994) models are only plotted 
to K — 17.8, but the star-count models and observations 
of Minezaki et al. (1998) are quite similar to the Reid et 
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, Star counts, this work 
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Figure 2. Rcdshift distribution predicted by the PLE model for 
all galaxies (solid) and E/SO galaxies only (dashed) to K' = 19.5. 



al. (1996) model for deep Galactic pole star-counts in the 
/-band, shifted by I - K ~ 1.25. 

We adopt the steeper Reid et al. (1996) /-band star- 
count model for the lower Galactic latitude of b = 45°, 
shifted by / — K — 1.25 and extrapolated at K < 14.5 with 
the Glazebrook et al. (1994) model for a b = 32° field. This 
is then normalized to our star-counts at 12.5 < K' < 17.0, 
with the best-fit normalization, corresponding to 1697 ±356 
deg _2 mag _1 at 17.0 < K' < 17.5. This model, shown on 
Figure 3, is then consistent with the observed counts within 
the statistical errors (reduced x 2 = 0.86). 



3 CLUSTERING AROUND THE RADIO 
GALAXIES 

3.1 Visual impression 

We first examine our data by eye for any obvious associa- 
tions of galaxies around the radio galaxies. Figure 4 shows 
maps of the distribution of K' < 19.5 galaxies on each of 
the five fields, with the radio galaxies indicated. The images 
provide more information than these maps as connecting 
filaments and interactions between galaxies may be visible. 
The visual impression is that 

(i) The radio galaxy 5C6.25 appears to lie on a long 
S-shaped filament, but as many of these galaxies are in fila- 
ments, this may not necessarily indicate an especially clus- 
tered environment. 

(ii) The radio galaxy 5C6.29 is not obviously in any sort 
of association (the four bright galaxies nearby are almost 
certainly at much lower redshift). 

(iii) 5C6.43 has a close companion of similar brightness 
(K' = 18.01), 6.4 arcsec away, with some signs of an inter- 
action, but it is not obviously part of any larger group. 

(iv) 5C6.75 lies within a more compact S-shaped struc- 



Figure 3. Differential star number counts in 0.5 mag intervals 
of K' magnitude, with field-to-field errors, and the adopted star- 
count model (dashed line). 



ture, some 1 arcmin across with at least 10 bright members 
and a number of faint galaxies in the same area. The region 
around the radio galaxy contains many interacting pairs and 
small groups of galaxies. The concentration of galaxies is ob- 
vious on Figure 3 and seems a good candidate for being a 
true cluster. 

(v) 6C0822 has two close companions of similar bright- 
ness and may be part of an interacting group of several 
galaxies, although with fewer members than the association 
around 5C6.75. 

Hence the radio galaxies appear to be in a wide range 
of environments, with evidence of clustering in some cases. 
To quantify this, we consider a model rich cluster placed at 
the redshifts of the radio galaxies. 

3.2 Cluster modelling 

One measure of the richness of a galaxy cluster is the 'Abell 
richness', Na, defined as the number of cluster members 
with apparent magnitudes < 2.0 magnitudes fainter than 
the third-ranked galaxy (e.g. Abell et al. 1989). Rich clusters 
are defined as those with Na > 30, with 30 < Na < 49 being 
Abell class 0, 50 < N A < 79 class 1 and 80 < N A < 129 
class 2. 

The surface density of galaxies in rich clusters, in terms 
of the projected distance from the centre r = 9<1a, typically 
follows a profile of the form 

Po 



p{r) = 



[1 + (r/Rc) 



where R c is the core radius, out to r„ 
zero. 



(1) 

where p(r) falls to 
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Figure 4. Positions of all A"' < 19.5 galaxies on the five Omega fields, each 5.75 X 5.75 arcmin. The short-dashed circles show a 1.5 
arcmin radius around each radio galaxy. The best-fit cluster centre on the 5C6.75 field is shown by an asterisk. 
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Following Kepner et al. (1998) we assume R c — 0.1 
h^ 1 Mpc (where h is Hubble's constant in units of 100 km 
s _1 Mpc _1 ) and a = 0.75, parameters in the middle of the 
observed ranges for both local and distant clusters (Girardi 
et al. 1995; Lubin and Postman 1996), and r max — 1.0h~ 
Mpc. With g = 0.05, a proper distance 1.0/i _1 Mpc at 0.7 < 
z < 0.8 corresponds to 3.38-3.57 arcmin, so the Omega field 
of view would contain almost the entirety of a cluster centred 
on a radio galaxy. 

Our model for field galaxies incorporates luminosity 
functions with a much steeper faint-end slope for blue than 
for red galaxies (see Roche and Eales 1998), but this may 
not be appropriate for rich clusters, which contain a much 
higher proportion of early-types and many low luminosity 
red galaxies. For example, Oemler et al. (1997) estimated 
that only 12 per cent of the galaxies in two high density 
Abell class 1 clusters at z ~ 0.4 were much bluer than pas- 
sively evolving ellipticals. To take this into account, when 
modelling the cluster we assumed the k+e-corrections of an 
evolving SO model for a large proportion (78 per cent) of 
the galaxies in the steeper luminosity functions correspond- 
ing to Sab and later types, to give a total red galaxy fraction 
~ 88 per cent. 

Kepner et al. (1998) parameterized the richness of 
model clusters in terms of the total luminosity (L c ; , in units 
of L*) within the assumed r max of 1/iT 1 Mpc, and estimated 
that this related to the Abell richness as Na — %L c i. We de- 
termine L c i for our model clusters by summing the luminos- 
ity of the cluster galaxies (to a faint limit K' — 22) relative 
to L* (hence L c i is not affected by luminosity evolution). 

For modelled clusters at the redshifts of the radio galax- 
ies, we estimate that ~ 0.7-0.8L c ; cluster members will be 
visible to the completeness limit of K' — 19.5, about half of 
these within 1.5 arcmin of the cluster centre. At z ~ 0.75, an 
L* elliptical [M K i = —24.85 for h=0.5) evolving as in our 
PLE model would have an apparent magnitude K' ~ 18.0 
and a typical brightest cluster galaxy {M K i — —26.49 for 
h = 0.5, Collins and Mann 1998) would have K' ~ 16.4, so 
any galaxies associated with the radio galaxies will almost 
certainly be fainter than K' — 16. We first estimate the 
cluster environment of the radio galaxies simply by count- 
ing the number of 16 < K' < 19.5 galaxies within 1.5 arcmin 
of the radiogalaxy positions and comparing with modelled 
clusters. 

3.3 Excess galaxies within 1.5 arcmin 

Table 2 gives the total number of 16 < K' < 19.5 galax- 
ies within 1.5 arcmin of each of the radio galaxy positions 
(including the radio galaxies themselves), N(< 1.5). The 
background galaxy density for each field is estimated by di- 
viding the number of galaxies > 1.5 arcmin from the radio 
galaxy, N(> 1.5) by the area of the part of the field > 1.5 ar- 
cmin from the radio galaxy, A(> 1.5). The excess E above 
the background of galaxies within 1.5 arcmin of the radio 
galaxy is first estimated by subtracting 

E« 1.5) = (N« 1.5) - ^ f 5) N(> 1.5) (2) 

However, if the radio galaxy does lie in a rich cluster, some 
cluster members - about half if the profile follows equation 
(1) - will lie > 1.5 arcmin from the cluster centre, causing 



Table 2. The number of 16 < K' < 19.5 galaxies within 1.5 
arcmin of the radio galaxy positions (within 1.0 arcmin in the case 
of 6C0822), the excess above the background density (equation 
2, with errors y/N ga i), the excess corrected for the presence of a 
cluster (equation 3), and the estimated total luminosity L c ; of a 
cluster centred on the radio galaxy. 



Galaxy 


N gal 


E(< 1.5) 


£/(< 1.5) corr 




5C6.25 


82 


11.5±9.1 


15.8 ± 12.5 


40 ±32 


5C6.29 


69 


8.2 ±8.3 


11.5 ± 11.6 


30 ± 31 


5C6.43 


62 


4.1 ± 7.9 


5.4 ±10.4 


15 ± 29 


5C6.75 


98 


33.5 ± 9.9 


45.0 ± 13.3 


125 ± 37 


6C0822 


36 


1.6 ±6.0 


2.0 ± 7.5 


8 ±28 



the background density to be overestimated. For a modelled 
cluster of richness L c i, we estimate a corrected background 
density by subtracting the number of model cluster galaxies 
> 1.5 arcmin from the cluster centre, C(> 1.5), from the ob- 
served N(> 1.5). This leads to a higher, corrected estimate 

E« 1.5) corr = N(< 1.5)- X> ls) ^ L5 )- C (> ^K 3 ) 

The cluster richness is then estimated as the value of L cl for 
which the number of galaxies within 1.5 arcmin of the centre 
of a modelled L c ; cluster, C(< 1.5), is equal to E(< 1.5) CO rr 
with the background correction from the same model clus- 
ter. For these fields, the background correction effectively 
increases the estimates of L c i by ~ 35 per cent. 6C0822 ac- 
tually gave a negative E(< 1.5) - it may coincide with a 
void in the distribution of nearer galaxies, so to reduce the 
effects of this we compute the excess for a smaller radius of 

1.0 arcmin. 

Of the 5 radio galaxies, only 5C6.75 shows a significant 
excess of nearby galaxies over the scales expected for a rich 
cluster. This excess would correspond to a cluster of Abell 
richness Na ~ 83 ± 25, i.e. Abell class 1 or 2. On the basis 
of y/N g ai statistics the cluster detection is 3.34<j, although 
this may be an overestimate due to the clustering of back- 
ground galaxies - see Sections 4.2 and 5.3 for estimates of 
the significance which take this into account. 

Of the other galaxies, 5C6.25 appeared to the eye to 
be in a large-scale association, but the estimated richness is 
probably too low for this to be described as even an Abell 
cluster, and 6C0822 appear to be in an association of a 
few galaxies, but the lack of any excess of galaxies over clus- 
ter scales of ~ 1 arcmin suggests this is a small isolated 
group rather than a true cluster. For the five radio galaxies, 
the mean L ct is 43.6 ± 21.1, corresponding to a mean Abell 
richness 29 ± 14 - this is discussed in Section 7.1. 

4 CLUSTER DETECTION 

4.1 Method 

We now apply a second method of searching for clusters, 
which makes use of the full model cluster profile of equation 
(1). This has the advantages of (i) providing an estimate of 
the cluster centroid position - assuming clusters are centred 
exactly on the radio galaxies is likely to underestimate the 
richness if there is an offset, (ii) providing an estimate of the 
richness which takes into account the full radial distribution 
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and should therefore give greater statistical power than a 
simple summation to a fixed radius, (iii) deriving errors from 
the real distribution of galaxies, rather than by assuming 
%JN statistics, and thus verifying that any cluster detections 
remain significant after taking into account the clustering of 
background galaxies. 

We search the regions near the radio galaxies for clus- 
ters using a cluster detection routine, which is essentially 
a simplified version of the Adaptive Matched Filter tech- 
nique described by Kepner et al. (1998). The profile of a 
cluster of chosen redshift z and richness L c i is modelled as 
in Section 3.2. For a chosen position (x,y), the number of 
16 < K' < 19.5 galaxies are counted in annular bins of 
A(8) = 10 arcsec centred on (x,y). The number of galax- 
ies in a model cluster, placed at (x,y), are counted in the 
same annular bins, with areas missed off the edges of the 
field or lost due to 'holes' taken into account. The model 
and observed counts are then compared by summing 



X (x,y,L c i,z) 



1 ^ 

i=l 



Aip 



(4) 



where d is the number of galaxies in bin i in the modelled 
cluster, Ni is the observed number of galaxies in bin i, Ai 
is the area of bin i (taking into account any area lost due 
to holes or the field edges), p the mean surface density of 
galaxies on the whole field and pback a corrected background 
density (the observed number of galaxies on the field minus 
the total number that would be in the modelled cluster, 
divided by the field area). 

The x 2 statistic is computed for a grid of (x,y) positions 
on each field, for a model cluster with z equal to the radio 
galaxy redshift and L c i — 125, the divide between Abell 
classes 1 and 2 and the approximate richness of the 5C6.75 
cluster. For each field u(x 2 ) is computed as the scatter in 
X 2 values between 4900 (x,y) positions, thus taking into ac- 
count both the non- independence of the bins and the greater 
than Poissonian variations in the galaxy surface density due 
to clustering. Candidate cluster positions are identified as 
minima in the x 2 map, where x 2 ls significantly, in terms of 
<t(x 2 ), lower than the mean for the whole field. 



4.2 Results 

The strongest detection of a cluster was centred at pixel co- 
ordinate (488,499) on the 5C6.75 field, where the x 2 for a 
L c i = 125 cluster profile was 3.06u below the mean for all 
positions on the field. This confirms that the cluster remains 
a 3<r detection after taking into account the increased varia- 
tions in the background density from the clustering of other 
galaxies on the field. There are no 3<r detections of clusters 
within 1.5 arcmin of any of the other four radio galaxies. 

The best-fit cluster position, shown on Figure 4, is off- 
set by 23 arcsec (~ O.llfe -1 Mpc) from the radio galaxy 
(478,442). However, running the cluster detection routine 
on 10 simulated L c i = 128 clusters, generated by distribut- 
ing galaxies within a probability distribution following the 
profile from Equation 1 on top of a randomly distributed 
background gives the RMS error in the best-fit position as 
~ 21 arcsec, so this offset is not significant. 

There are 99 galaxies with 16 < K' < 19.5 within 1.5 
arcmin of the fitted cluster centre, an excess of 34.4 ± 9.9 



1 i , I 1 



100 150 200 

Distance from centre (arcsec) 



Figure 5. Radial density profile of the 5C6.75 cluster, with y/N 
error bars, compared to (dotted line) a model cluster profile 
(equation 1) with L c i = 128 and R c = 0.1 h^ 1 Mpc. 



above the background or 46.6 ± 13.4 above the corrected 
background, giving L c i = 128 ± 37 (hence Abell richness 
Na — 85 ± 25), almost identical to the estimate when the 
radio galaxy is assumed to be the centre. The richness can 
also be estimated using the cluster detection routine, placing 
a model cluster at the best-fit position and finding the model 
L c i which mimimizes x' ' ■ This gives L c i = 136 ± 26 (error 
from the simulation), consistent with the first estimate. 

Figure 5 compares the observed cluster profile with the 
Lci = 128, R c = 0.1 h' 1 Mpc model. Allowing both L c i 
and R c to vary gives a best-fitting core radius as R c — 
0.12 ± 0.03fe _1 Mpc. The good agreement of the observed 
cluster profile and its best-fit R c with the model profile con- 
firms that the distribution of galaxies near 5C6.75 is at least 
consistent with an Abell 1/2 cluster, similar in profile to 
those in the local Universe, at the radio galaxy redshift. 



5 THE ANGULAR CORRELATION 
FUNCTION 

5.1 Calculating lo(8) 

We investigate the clustering of all the galaxies on these five 
fields by calculating the angular correlation function, u>(9). 
Our data may not give an entirely unbiased estimate of the 
field galaxy u>(6), due to the presence of known radio galax- 
ies with strong clustering around at least one. However, the 
great majority of detected galaxies will not be associated 
with the radio sources, so the resulting bias may be very 
small. For now, we analyse these images as normal field sam- 
ples, with no distinction between radio and other galaxies, 
but investigate the effect of the 5C6.75 cluster at the end of 
Section 5.2. 

For each field, u>{6) is calculated for all detections 
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classed as galaxies to a series of faint magnitude limits. On 
a field with N g galaxies brighter than the chosen limit, there 
will be 4\N g (N g — 1) possible galaxy-galaxy pairs, which 
counted in bins of separation of width A(log 9) — 0.2, giv- 
ing a function N gg (9i). For each field N r — 10000 points are 
placed at random over the area covered by real data, i.e. 
avoiding any holes, and the separations of the N g N r galaxy- 
random pairs, taking the real galaxies as the centres, are 
similarly counted, giving N gr (9i). In addition, the separa- 
tions of the ^N r (N r — 1) random-random pairs are counted 
in bins to give N rr (6i). 

Defining DD — N gg (9i), and DR and RR as the galaxy- 
random and random-random counts normalized to have the 
same summation over all 9 as DR, i.e. DR = ^ N ^ N ^ N gr (9i) 

and RR = (jv° _ i) N rr (9i), we calculate uj(9) for each 0, 
bin using the Landy and Szalay (1993) estimator, 



DD - 2DR + RR 
RR 



(5) 



The u>(6i) of the five fields are averaged at each 9i point 
to give a mean u>(9), shown on Figure 6 for magnitude lim- 
its K' = 18.5-20.0, with error bars from the field-to-field 
scatter. 

If the real galaxy u>(9) is of the form A9~ s , the observed 
u(9) will follow the form u(9) = A(9~ s -C), with amplitude 
A (defined here at a one-degree separation), and a negative 
offset AC (the integral constraint) resulting from the re- 
stricted area of the observation, which can be estimated by 
doubly integrating an assumed true u>(6) over the field area 

n, 



AC = 



Q 2 



Lu(9)dn 1 dn 2 



(6) 



Using the random-random correlation, this calculation can 
be done numerically - 

£iV rr (0)0- a 



C 



(7) 



T,N rr (9) 

Assuming S — 0.8, in agreement with most observations in- 
cluding those in the if-band at brighter limits (Baugh et 
al. 1996), C = 14.84 for these fields. The uj(9) amplitude A 
is then estimated by fitting A{0 — 14.84) to the mean 
to(9) at separations 2 < 9 < 200 arcsec. The error on A is 
estimated by fitting the same function to the u>(9) of the 
five individual fields and determining the scatter between 
the individual field amplitudes. 

5.2 uj(9) results 

Figure 6 shows the observed uj(9) with best- fitting functions 
of the form uj{9) = A(9' - 8 - 14.84), with the amplitudes 
A in Table 3. At magnitude limits K' = 18.5-20.0, there is 
a significant (~ 4er) detection of galaxy clustering, with the 
individual A(log 9) = 0.2 bins showing a significant (~ 3u) 
and positive signal out to separations log 9 = —2.25 = 20 
arcsec. The result for K' — 20 may be less reliable due to 
incompleteness at K' > 19.5 (Section 2.3), whereas at limits 
of K' — 18.0 and brighter, there are an insufficient number 
of galaxies for a significant detection of clustering. 

The first two plotted points, corresponding to 1.26 < 
9 < 3.17 arcsec, lie above the fitted power-laws, indicating 
an excess of close pairs of galaxies relative to the clustering 



Table 3. Observed w(6) amplitudes (A), in units of 10~ 4 at 
1°, the number of galaxies (N ga i), the estimated fraction of con- 
taminating stars f s and the ui(6) amplitudes corrected for star- 
contamination A corr , at a series of K' magnitude faint limits. 



K' limit 


A 


N gal 


/. 


A C orr 


18.5 


12.77 ± 1.84 


807 


0.180 ±0.038 


21.94 ±3.25 


19.0 


11.28 ±2.98 


1117 


0.164 ±0.034 


16.14 ±4.47 


19.5 


13.25 ±2.27 


1626 


0.150 ±0.032 


18.34 ±3.43 


20.0 


11.20 ±2.19 


2110 


0.148 ±0.031 


15.43 ± 3.22 



seen at larger scales. In Paper I, the large number of 2-3 arc- 
sec pairs on the Redeye camera fields was attributed to the 
effects of galaxy mergers, but it was also suggested that u>{9) 
from Jf-band surveys might be steeper than 9~°' s , due to 
a particularly steep two-point correlation function, £(r), for 
giant ellipticals. Loveday et al. (1995) measured a £(r) slope 
7 = 1.87 ± 0.07 for E/S0 galaxies and Guzzo et al. (1997) a 
more extreme 7 = 2.05 ± 0.09 for higher luminosity E/SOs. 
The latter estimate corresponds to co(9) oc 9~ 105 , which we 
assume as an extreme upper limit to the possible range of 
slopes, and fit the observed ui{9) with i A m {0~ 1 ' 05 — C m )\ 
where C m = 38.19 for a -105 power-law (from equation 7). 

Figure 6 shows this steeper power-law with the best- 
fitting amplitudes (A m = 3.93 ± 0.39 x 10~ 4 , 2.88 ± 0.62 x 
10~ 4 , 3.132 ± 0.49 x 10~ 4 , 2.49 ± 0.46 x 10~ 4 at limits 
K' = 18.5, 19.0, 19.5 and 20.0 respectively). Changing from 
9~ 8 to #~ 105 does not get even half-way to fitting the 
small-scale excess (except perhaps at the K pr " lme — 20 limit 
where close-pair detection is likely to suffer incompleteness) , 
suggesting that the dominant contribution to this is from in- 
teracting galaxies rather than steeper E/S0 clustering. The 
excess of close pairs is investigated in Section 6 and discussed 
in Section 7.3. 

A fraction of (randomly distributed) stars f„ within the 
galaxy sample will reduce the observed u>(6) at all angles 
by a factor of (1 — f s )~ 2 - For each magnitude limit, we esti- 
mate, using the star-count model of Section 2.3, the fraction 
f s of faint (K > 17) stars contaminating the sample of N ga i 
objects. Table 2 gives f s and ui(9) amplitudes A corr cor- 
rected for star contamination by multiplying the fitted A by 
(l~fs)- 2 , With errors combining in quadrature the co(9) er- 
rors with the correction error (derived from the 21 per cent 
uncertainty in the star-count normalization). 

It might be thought that the cluster on the 5C6.75 field 
would upwardly bias our estimate of the field galaxy w{9), 
but repeating the analysis with galaxies within 1.5 arcmin 
of the cluster centre excluded produces little change in the 
results (in units of 10~ 4 , the uncorrected u>(8) amplitudes 
become 10.10±2.55, 10.09±2.53, 15.13±2.85 and 11.21±2.73 
at K' = 18.5,19.0,19.5 and 20.0 respectively). Although 
the £(r) amplitude in the core of Abell 1/2 clusters may be 
several times higher than for field galaxies (see e.g. Hill and 
Lilly 1991), the excess galaxies within 1.5 arcmin of 5C6.75 
amount to only 13 per cent of all galaxies on the CCD frame, 
which in turn is only one-fifth of our dataset, so the effect 
of their clustering on u)(0) is greatly diluted. 
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Figure 6. Observed u>{6) for galaxies brighter than K' = 18.5, 19.0, 19.5 and 20.0, as log-linear (left) and log-log (right) plots. The 
dashed lines show the best-fit uj(8) = A(8~ ihR — 1.84 functions and the dotted lines the best-fit functions for a maximum slope, of the 
form A m (0- lm - C). 
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Figure 7. The ui(0) amplitude for galaxies on our five fields, 
with and without corrections for star-contamination, against K' 
magnitude limit, compared with u)(Q) amplitudes from the K- 
band surveys of Baugh et al. (1996), Carlbcrg et al. (1997) and 
Paper I. The solid line shows a PLE model (Roche and Ealcs 
1998) with ro = 8.4h Mpc for early-type galaxies and stable 
clustering (e = 0), the dashed line the same PLE model with 
e = 1.2 evolution of the clustering, the dot-dash line the same 
model with increased clustering at z > 1.5 (sec Section 7.2), and 
the dotted line the e = PLE model with the ro from Roche et 
al. (1996), r = 5.9/1- 1 Mpc for E/SO galaxies. 



Figure 7 shows the corrected and uncorrected w{ff) am- 
plitudes as a function of K' limit, from this survey and oth- 
ers in the K band, and models described in Section 5.4. 
Note that the Paper I and Carlsberg et al. (1997) results are 
plotted without a correction for star-contamination (Baugh 
et al. 1996 are at sufficiently bright magnitudes that star- 
galaxy separation should be reliable). At K = 19.5-20.0, our 
uncorrected amplitudes are consistent with Paper I. 



5.3 Effect on detection of rich clusters 

Having directly measured the galaxy clustering on these 
fields we can estimate its effect on the detection of rich clus- 
ters around the radio galaxies. Considering cells of area Q, 
with a mean background density of (N) galaxies in each, 
clustering will increase the variance ^2 of the counts in these 
cells as 



= (N)+ = (N)' 



(8) 



AO and a circular area of radius 1.5 

2 



which for ui(9) = 

arcmin gives /12 = 25. 6 A(Ny . The background density 
for the 5C6.75 field was estimated as (N) = 53 and the 
uncorrected (as the background density will include con- 
taminating stars) A = 1.325 x 10~ 3 , giving a background 
a = ^/(53 + 95.3) = 12.2. The significance of the excess of 
galaxies within 1.5 arcmin of 5C6.75 is then reduced from 



3.34,7 to f|| = 2.75a. This is slightly lower than the 3.04a 
significance from the cluster detection routine, which is pre- 
sumably more statistically powerful as it takes into account 
the full cluster profile. 

We conclude that the effects of the clustering of unas- 
sociated galaxies on our investigation of radio galaxy cluster 
environments are non-negligible but relatively small, a ~ 20 
per cent increase in statistical errors, and that the detection 
of the 5C6.75 cluster remains close to 3a. 



5.4 Comparison with models 

The modelling of co(8) amplitudes is disussed in detail else- 
where (e.g. Roche et al. 1996; Roche and Eales 1998) and 
will be described only briefly here. Essentially, the 3D two- 
point correlation function for the galaxies is parameterized 



ar,z) = (r/r )- 1 (l + z) 



"(3+0 



(9) 



where ro normalizes the strength of clustering at z = 0, 7 
is the slope and e the clustering evolution relative to the 
e = stable clustering model. This is integrated over a 
galaxy redshift distribution, N(z), using Limber's formula 
(e.g. Phillipps et al. 1978). Here N(z) is given by the PLE 
model of Section 2.3. 

Roche et al. (1996) assumed 7 = 1.8 for all galaxy types, 
with ro = 5.9 h^ 1 Mpc for E/S0 galaxies and r = 4.4 h^ 1 
Mpc for later types, with an additional luminosity weight- 
ing so that galaxies with 2 = blue-band absolute magni- 
tudes M B > -20.5 (for H = 50 km s _1 Mpc _1 ) are half as 
clustered as those of higher luminosity (see Loveday et al. 
1995). With stable clustering and L* evolution, this model 
predicted a u>(6) scaling consistent with observations in the 
blue-band, but appeared to underpredict the lo(8) ampli- 
tudes from K-b&nd surveys at K ~ 19.5-21.5 (Paper I). The 
results from the deep if-band surveys, in which early-type 
galaxies will be much more prominent, were better fitted by 
a model with even stronger clustering for early-type galax- 
ies. 

The u){9) scaling from a large i?-band survey (Roche 
and Eales 1998) appeared well-fitted by a PLE model with 
7 = 1.8, r = 8.4 h^ 1 Mpc for E/S0 galaxies and r = 4.2 
h^ 1 Mpc for spirals and irregulars, again with the same lu- 
minosity weighting, and e ~ 0, whereas clustering evolution 
of e ~ 1.2 underpredicted the lo{9) amplitudes. On Figure 
7, the same PLE model with e = but the ro from Roche 
ct al. (1996) underpredicts the co(8) amplitudes from this 
survey and is rejected after star-contamination is taken into 
account. With the ro from Roche and Eales (1998), the PLE 
model with e = gives higher u(0) amplitudes more consis- 
tent with the observations. The same model with clustering 
evolution of e = 1.2 greatly underpredicts the u>(9) of deep 
K-band surveys, and at K' — 19.5 is rejected by > 3a even 
before the correction for star-contamination. 

At our fainter limits of K' = 19.5-20.0, there may be 
some evidence that the w(0) amplitude falls less steeply than 
the e — model, and the amplitude from the small area 
studied by Carlberg et al. (1997) K' — 21.5 suggests that 
the u>(6) scaling actually levels out. This might indicate an 
increase in the clustering of red galaxies at the highest red- 
shifts (z > 1.5), perhaps related to the enhanced clustering 
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reported for Hubble Deep Field galaxies at z > 2.4 (Maglioc- 
chetti and Maddox 1998), but more data is needed to con- 
firm this. If we assume the levelling out of the u>(9) scaling 
is genuine, it can be fitted by the e = PLE model with 
an increase in clustering to ro = 16.15/i _1 Mpc (see Section 
7.2) for all z > 1.5 galaxies, as shown on Figure 6. 

In summary, the uj (9) amplitudes from this survey are 
consistent with previous K-b&nd observations, favouring 
both strong (r () ~ 8.4 h^ 1 Mpc) clustering for E/SO galaxies 
and clustering remaining stable (e ~ 0) to z ~ 1.5; if any- 
thing the clustering is stronger than an e = model at the 
highest redshifts. We discuss this further in Section 7.2. 



6 CLOSE PAIRS OF GALAXIES 
6.1 Method 

The lj(9) from this data exceeded the fitted power- law at 
small separations of 1.26 < 6 < 3.17 (Figure 6), indicating 
that there is an excess of close pairs of galaxies. As in paper 
I, we quantify this using a method described by Woods, 
Fahlman and Richer (1995). A probability P of occurring by 
chance (in a random distribution of galaxies) is estimated 
for each galaxy-galaxy pair, as 



P= ff exp(— Tvpa 2 )da 



(10) 



where p is the surface density of galaxies brighter in apparent 
magnitude than the fainter galaxy of the pair, 9 is the pair 
separation and /3 an angular separation cut-off below which 
individual objects cannot be resolved. 

We take (5 = 1 arcsec and find the area around each 
galaxy by counting randomly distributed points in annular 
bins rather than assuming irr , thus taking into account field 
edges and holes. For all 1629 objects classed as galaxies to 
the estimated completeness limit K' = 19.5, the pairs with 
P < 0.05 are counted in A(9) = 0.5 arcsec bins of separa- 
tion. We then perform an identical analysis on 25 random- 
ized datasets in which the same number of 'galaxies', with 
the same magnitudes, are redistributed randomly over the 
field areas. The pairs counts from the randomized dataset 
and averaged, and their scatter used to derive error bars for 
a single dataset. 



6.2 Results 

Figure 8 shows a histogram for the P < 0.05 pair counts 
compared to the random expectation. Over the 1.5 < 9 < 
5.0 arcsec range, there are 157 pairs with P < 0.05 com- 
pared to 63.48 ± 10.57 expected for a random distribution 
of galaxies. The excess above random, assuming the error to 
scale as y/N , consists of 93.52 ± 16.62 pairs. 

One one field, every pair was examined by eye, and only 
2 out of 38 appeared spurious (caused by irregular outer 
regions of bright galaxies), while the other 36 appeared to be 
genuine galaxy pairs, in many cases with signs of interaction. 
There are 15 1.5 < 9 < 5.0 arcsec pairs within 1.5 arcmin 
of the cluster radio galaxy 5C6.75. This can be compared to 
3.16 expected for a random distribution, increasing to 7.8 or 
10.8 if, respectively, the pairs/randoms and pairs/galaxies 
ratios for this area are the same as for the whole dataset. 
Hence our discovered cluster may contribute a few of the 
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Figure 8. Histograms of the observed numbers of galaxy-galaxy 
pairs (to K' = 19.5) with probability P < 0.05. The dotted 
line shows the pair count from simulated random distributions 
of galaxies, with error bars from the simulations. 



excess pairs, but increases the estimate of the fraction of 
galaxies in pairs by only ~ ^ 15 ^ 2 4 ^ = 2.7 ± 1.3 per cent. 

No 9 < 1 arcsec pairs are counted, as with the FWHM ~ 
1.13 arcsec resolution they would be merged into single de- 
tections. The observed number of 1.0 < 9 < 1.5 arcsec pairs 
is approximately the random expectation, suggesting (by 
comparison with the number of 1.5 < 9 < 5.0 arcsec pairs) 
that only ~ 40 per cent of the true number are resolved. The 
number of close pairs missed due to image merging is esti- 
mated by assuming that all 9 > 1.5 arcsec pairs are resolved 
and that the small-scale lo(9) is still a #~ 0,8 power-law, which 
gives the true number of pairs above the random expecta- 
tion at 9 < 1.5 arcsec as 0.3068 the number at 1.5 < 9 < 5.0 
arcsec. 

To estimate the number of physically merging or inter- 
acting galaxies, we must subtract from the observed num- 
ber of pairs, the number expected from normal galaxy clus- 
tering, which will be higher than the random expectation. 
Clustering with a 9~ o a power-law of amplitude a at 1 arc- 
sec increases the ratio of 9 < 5 arcsec pairs to the random 
expectation by 



f 5 2n9a9- 8 d9 _ 2to fl 1 - 2 
~ 25:; 1.2 



0.46a 



(11) 



77(5)2 

Our measured co(9) amplitude at the K' = 19.5 limit cor- 
responds to a = 0.93 ± 0.16 without correction for star- 
contamination or ex = 1.28 ± 0.24 with correction, and the 
random expectation for 9 < 5 arcsec pairs is 79.28. Hence, if 
we first neglect the star-contamination, the excess of 9 < 5 
arcsec pairs above the expectation from the larger scale lo(9) 
amounts to 1.3086 x (93.52 ± 16.62) - 0.46 x (0.93 ± 0.16) x 
79.28 = 88.46 ± 22.52. As each pair consists of two galax- 
ies, the fraction of the 1629 galaxies in merging/interacting 
pairs is then (2 x 88.46 ± 22.52)/1629 = 10.9 ± 2.8 per cent. 
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Taking into account the star contamination, the excess 
of 9 < 5 arcsec pairs is reduced to 1.3086 x (93.52 ± 16.62) — 
0.46 x (1.28 ±0.23) x 79.28 = 75.70 ±23.31, but the number 
of real galaxies is also reduced, to 1382, leading to a very 
similar estimate of the merging/interacting fraction as (2 x 
75.70 ± 23.31)/1382 = 11.0 ± 3.4 per cent. This estimate is 
discussed further in Section 7.3. 



7 DISCUSSION 

7.1 The radio galaxy clustering environment 

At low redshifts, the most radioluminous (FRII) radio galax- 
ies tend not to be found in galaxy clusters, whereas the 
much lower luminosity FRI sources often are. This scenario 
changes at relatively moderate redshifts - Hill and Lilly 
(1991) found that 43 radio galaxies at 0.35 < z < 0.55, with 
a very wide range of radio luminosities, were distributed in 
similar numbers over field, Abell and Abell 1 cluster en- 
vironments. The mean environment of the 43, quantified in 
terms of B gg , the normalization of £(r) at r = 0.5/t -1 Mpc, 
was estimated as B gg — 291 ± 45, approximately an Abell 
cluster, with no strong dependence on radio luminosity - for 
the most luminous (3C) subsample the mean environment 
was B gg = 342 ± 96. It was concluded that the mean B gg of 
luminous (FRII) radio galaxies (including those of with the 
radio luminosities of our 7C sample) increased by a factor 
~ 2.5 from z ~ to z ~ 0.5. We aim to determine whether 
there is a further increase at z > 0.5. 

Previously, Bower and Smail (1997) detected by grav- 
itational lensing a cluster of apprximately Abell 2 richness 
around a radio-loud QSO (3C336) at z = 0.927, but found 
no significant lensing signal for 7 other 3C radio galaxies at 
0.87 < z < 1.05, indicating a mean environment no richer 
than an Abell 1 cluster. In Paper I, we were able only to set 
a similar upper limit (B gg < 756 at 2a) on the strength of 
clustering around radio galaxies at z mean ~ 1.1, but here, 
due to the more moderate (0.7 < z < 0.8) radio galaxy red- 
shifts and the larger field size, we are more successful. We 
detect at ~ 3a significance a cluster of estimated Abell rich- 
ness Na = 85 ± 25 (class 1 or 2), approximately centred on 
the radio galaxy 5C6.75, with a profile consistent with that 
of a typical present-day rich cluster (R c ~ 0.1 h^ 1 Mpc, 
a ~ 0.75 in equation 1) at the radiogalaxy redshift. Of the 
other four radio galaxies, two appeared to be in groups or 
structures of lower richness than Abell clusters, and two in 
field environments. 

We estimated the mean environment of our sample of 
0.7 < z < 0.8 radio galaxies as Abell richness Na = 29 ± 14. 
Using large samples of clusters, Lubin and Postman (1997) 
find N A ~2.5-2.7 JV , where N is the core (< 0.25 h' 1 
Mpc) richness, and Hill and Lilly (1991) fit B gg ~ 30iVo. 
Combining these two relations, Na = 29 ± 14 corresponds 
to B gg = 335± 162, consistent with the Hill and Lilly (1991) 
estimate for radio galaxies at 0.35 < z < 0.55. Hence, al- 
though a larger sample will be required to confirm this, our 
results together with Paper I and Bower and Smail (1997), 
suggest there is little further evolution in the distribution of 
luminous radiogalaxy environments from z = 0.5 to z = 1. 

Wan and Daly (1996) explain the increase in the mean 
richness of FRII radio galaxy environments between z — 



and z — 0.5 as the result of an increase with time in the 
typical intra-cluster medium pressure; once the pressure in 
a cluster reaches a critical threshold FRII activity is sur- 
pressed and only the much less radioluminous FRI outbursts 
will occur, so FRIIs tend not to be found in clusters at 
lower redshifts. We infer from this that, firstly, if FRI radio 
galaxies are not surpressed by intra-cluster pressure, their 
present-day distribution - a wide range of field and cluster 
environments with the mean richness being about that of 
an Abell class cluster (Prestage and Peacock 1988; Hill 
and Lilly 1991; Zirbel 1997) - should reflect the range of 
environments of massive black holes capable of radio out- 
bursts. Secondly, if at some high redshift, all clusters have 
an intra-cluster pressure below the threshold for FRII sur- 
pression, the distribution of FRII environments should then 
be similar to that of the low redshift FRIs. 

Hence, on the basis of the Wan and Daly (1996) model, 
the rising trend in the mean richness of FRIIs environments 
with redshift should approach - but not exceed - that of 
an Abell cluster at high redshifts, and, as the Abell 
richness has already been reached at z ~ 0.5, there would 
be little further change at z > 0.5. This is consistent with 
our observations. 

The removal of a strong cluster-environment influence 
on the activity of radio galaxies beyond z ~0.5-0.6 might 
also help to explain the tighter correlation of radio luminos- 
ity to host galaxy mass for z > 0.6 radio galaxies, compared 
to those at lower redshifts (Eales et al. 1997; Roche, Eales 
and Rawlings 1998). 

7.2 The galaxy ui(9) amplitude 

At magnitude limits K' — 19-20, we detect galaxy cluster- 
ing at ~ 4a significance. In Paper I we interpreted the rela- 
tively high uj{9) amplitude at K ~ 20 as evidence that the 
red early-type galaxies were much more clustered than spi- 
rals, even at z ~ 1, and followed a stable clustering (e ~ 0) 
model. This paper's results support these conclusions at bet- 
ter statistical significance. The if-band lo(9) scaling is rea- 
sonably consistent with a PLE model with e = and ro of 
4.2k' 1 Mpc for spirals, 8.4/2T 1 Mpc for E/S0 galaxies. The 
high lo(9) amplitude and its slow decline on going faintward 
strongly favour e ~ and disfavour clustering evolution of 
e ~ 1.2 by ~ 3a. 

The larger R < 23.5 survey of Roche and Eales (1998) 
found uj{9) to be consistent with the same PLE model as 
used here, with stable clustering, and Postman et al. (1998) 
found the ui(9) in even larger 7 < 23 survey best-fitted by an 
essentially similar model with N(z) derived from the CFRS, 
e = -0.20 ± 0.17 and a mean r = 5.6 ± 0.23ft _1 Mpc. As 
most of the higher redshift galaxies in a 2.1/im survey will be 
red E/SOs (Figure 3), whereas the R and 7- band surveys will 
be dominated by later-type galaxies, the fact that surveys 
of similar depth in R, I and K' all favour stable clustering 
suggests e ~ for the separate E/S0 and spiral populations. 
The former would then be ~ 0.5 dex more clustered at all 
redshift (to at least z ~ 1.5), as suggested by the Neuschae- 
fer et al. (1997) u>(0) analysis where bulge-profile galaxies 
on deep HST images remain ~ 0.5 dex more clustered than 
disk galaxies from I ~ 18 to I ~ 24. 

CDM models with Q = 1 predict e ~ 1—1.5 for the mass 
distribution over the redshift range of our survey (Col/'in, 
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Carlberg and Couchmann 1997; Moscardini et al. 1998), 
whereas the inferred e ~ is more consistent with low Q, 
cosmologies. In the models of Moscardini et al. (1998), e ~ 
for galaxies from z ~ to z ~ 1.5 would only be consistent 
with Q = 1 if there is a rapid increase of the biasing of 
galaxies relative to the mass distribution over this redshift 
range, but this possibility is found to be excluded by the 
low u){9) amplitude of the Hubble Deep Field and Canada- 
France Redshift Survey galaxies. 

Enhanced bias is only seen on the /-limited HDF at 
z > 2.4 (Magliocchetti and Maddox 1998), which is proba- 
bly beyond the redshifts reached by our survey (Figure 2). 
However, if the Lyman Break Objects subsequently become 
passively evolving red galaxies at z ~ 1.5-2, with very red 
colours of I — K' ~ 5, they would then be much more promi- 
nent in a deep if-band sample than at I and shorter wave- 
lengths, and the influence of their strong clustering might 
extend to lower redshifts. We consider the possibility that 
this effect is first seen at our fainter limits of K' = 19.5-20, 
where the ui(9) amplitudes exceed the e = PLE model, 
becoming larger at the K — 21.5 limit of Carlberg et al. 

(1997) . 

The clustering measured by Magliocchetti and Maddox 

(1998) for HDF galaxies with photometric redshifts 2.4 < 
z < 2.8 corresponds, for a low Q, to a comoving ro = 6.875 ± 
2.7h~ 1 Mpc, which at z = 2.6 is equivalent to a stable model 
with ro = 16.15 ± 6.4ft -1 Mpc. The dot-dash line on Figure 
7 shows the e = PLE model in which the E/SO and spiral 
ro are both increased to this high value above a redshift Zi. 
To fit the high uj(6) amplitudes at K' > 19.5 requires Zi ~ 
1.5 for /^-limited surveys, rather than the Zi ~ 2.4 directly 
observed in the /-limited HDF. Consistency of this scenario 
with the much lower u)(6) amplitudes seen in the deepest A ~ 
0.4-0.9/xm surveys would require there to be a very strong 
correlation between galaxy colour and cluster environment 
at z ~ 1.5-3. This might result from an extended epoch 
of galaxy formation in which galaxies formed significantly 
earlier, and evolved with shorter star-formation timescales, 
in denser regions of the Universe. 

7.3 The evolution of the galaxy merger rate 

Infante et al. (1996) reported that, for a sample of R < 21.5 
galaxies, the uj(6) at 9 < 6 arcsec exceeded the inwards 
extrapolation of the uj(9) at large separations by about a 
factor of 1.8. This excess of close pairs was interpreted as 
consisting of physically interacting or merging galaxies, with 
the numbers suggesting some increase with redshift in the 
merger/interaction rate. In Paper I, there were more than 
twice as many 2 < 8 < 3 arcsec separation galaxy pairs 
as expected by chance for the section of the data observed 
with the Redeye camera, but the excess of pairs above the 
higher expectation from the fitted u>(9) was of marginal sig- 
nificance, and no close pair excess was seen on the other 
half of the data observed with the Magic camera. In the 
much larger _R-band survey of Roche and Eales (1998), one 
of the observed fields showed a significant excess in ui(6) at 
2 < 8 < 5 relative to the amplitude at 9 > 5 arcsec, again 
suggesting merger-rate evolution, but the second field stud- 
ied did not show a significant difference in the small and 
large scale ui(6) amplitudes. 

On the five fields studied here, the observed number of 



close (1.5-5.0 arcsec separation) pairs of galaxies exceeds, 
at ~ 3.3<r significance, the expectation from the fitted ui(9). 
Similarly, the to(9) at 6 < 3.2 arcsec shows an excess over the 
fitted power-law (Figure 6), even if u>(9) even for a steeper 



power- law of to(9) 



corresponding to the steepest 



f (r) observed for any type of galaxy (Guzzo et al. 1997). 

After adding an estimate of the number of unresolved 
9 < 1.5 arcsec pairs, we estimate that 11.0 ± 3.4 per cent 
of K' < 19.5 galaxies belong to 9 < 5 arcsec pairs in excess 
of the fitted u>(9), if a #~ 0,8 power-law is assumed. This is 
consistent with the estimated pair fractions on the Redeye 
fields of Paper I at K < 20 and on field 'e' of Roche and Eales 
(1998) at R < 23. The clear detection of a close pair excess 
in this relatively high resolution (FWHM ~ 1.13 arcsec) 
data suggests that the non-detection on the Magic fields of 
Paper I and field 'f of Roche and Eales (1998) was simply 
the result of the poorer resolution (FWHM ~ 1.8 arcsec) of 
these datasets. 

The fraction of galaxies in close pairs can be interpreted 
in terms of an evolving merger/interaction rate, parameter- 
ized here as R m (z) oc R m o(l + z) m . Following Roche and 
Eales (1998), we assume f pa ir(z) oc R m (z) (on the basis that 
the merger timescale is much shorter than the Hubble time) 
and normalize f pa ir at z — using the Carlberg et al. (1994) 
estimate that 4.6 per cent of local galaxies are in pairs of 
projected separation < 19 h^ 1 kpc. All such pairs will have 
an angular separation 9 < 5 arcsec at angular diameter dis- 
tances (La > 784 h^ 1 Mpc, i.e. at z > 0.455 in our chosen 
cosmology. At lower redshifts, assuming lj(9) oc #~ 0,8 , the 
fraction of these close-pair galaxies with 9 < 5 arcsec will be 

/ d A (z) -.1.2 



We then model f pa ir by summing 



f P air(z < 0.455) + f palr (z > 0.455) 



(12) 



where 



/„<,«■(« < 0.455) = 0.046- 



fparr(z > 0.455) = 0.046 



j;r 5 N(z) dz 

/ 6 455 ^(*)( 1 + *)"' d * 
/o 6 465 



over the PLE model N(z) for K' < 19.5 (Figure 3). This 
predicts f pa i r = 4.3 per cent for no evolution (m = 0), 
8.6 per cent for m = 1 and 18.2 per cent for m = 2, with 
the observed f pair — 11.0 ± 3.4 per cent corresponding to 
m = 1.33lS;£. 

This merger rate evolution can be compared with the 
m = 2.2 ± 0.5 estimated by Infante et al. (1996) from R < 
21.5 galaxies, the m = 2.01± \H of Roche and Eales (1998) 
for the R — 21.5-23.5 range of limits, and the m = 1.2 ± 
0.4 of Neuschaefer et al. (1997) from deep (I < 25) HST 
data. Taken in combination, these estimates may hint at a 
reduction in m with increasing depth. 

Carlberg et al. (1994) predict that m will be sensitive 
to f2 but will also be affected by any reduction in the mean 
mass of galaxies at higher redshift (inevitable if merging 
is occurring), which would cause the merger rate evolution 
to depart from a simple R m (z) oc (1 + z) m form in that 
deeper surveys would give lower estimates of m. For exam- 
ple, an Q — 0.18 Universe is predicted to give m = 2.2 at 
low redshift, falling to m = 1.2 at the redshift where the 
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mean galaxy mass is reduced by 32 per cent. Our results 
and the others discussed here would be consistent with this 
scenario, but inconsistent with the more rapid merger-rate 
evolution (m = 3.2-4.5) expected for Q. = 1. To accurately 
determine the form of R m (z), large and deep surveys with 
multiple passbands (ideally both the near-infra-red and op- 
tical) to provide photometric redshift estimates will be re- 
quired. Studies of close pair statistics appear to be critically 
dependent on data resolution, so any attempt to determine 
Rm(z) should use either space-based observations or only 
the sharpest of ground-based images. 
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